Intra-cellular traffic: bio-molecular motors on filamentary tracks* 
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Molecular motors are macromolecular complexes which use some form of input energy to perform 
mechanical work. The filamentary tracks, on which these motors move, are made of either proteins 
(e.g., microtubules) or nucleic acids (DNA or RNA). Often, many such motors move simultaneously 
on the same track and their collective properties have superficial similarities with vehicular traffic 
on highways. The models we have developed provide "unified" description: in the low-density limit, 
a model captures the transport properties of a single motor while, at higher densities the same 
model accounts for the collective spatio-temporal organization of interacting motors. By drawing 
analogy with vehicular trafflc, we have introduced novel quantities for characterizing the nature of 
the spatio-temporal organization of molecular motors on their tracks. We show how the traffic-like 
intracellular collective phenomena depend on the mechano-chemistry of the corresponding individual 
motors. 

PACS numbers: 89.20.-a;89.75.-k 



I. INTRODUCTION 

Motility is the hallmark of life. A distinguishing fea- 
ture of molecular motor transport in eukaryotic cells is 
that the motor proteins move on filamentary "tracks" 
[11, S, S 0] ■ The tracks for motor proteins are made of 
either proteins or nucleic acids. Not all motor proteins 
carry molecular cargo. A common feature of all these 
motors is that these perform mechanical work by utiliz- 
ing some other form of input energy and hence the name 
"motor" . All the molecular motors we consider in this 
paper directly convert chemical energy into mechanical 
work. 

During several biological processes many motors move 
simultaneously on the same track. The collective move- 
ment of the motors under such cicumstances strongly 
resemble vehicular traffic flow ^l- In this paper we 
present few models of molecular motor traffic which, un- 
like earlier works of other groups, also capture the essen- 
tial steps of the mechano-chemistry of individual motors. 

Our models of molecular motor traffic are biologically 
motivated extensions of the totally asymmetric simple 
exclusion process (TASEP) [7]. A TASEP is defined in 
terms of "rules" for updating the states of the system: 
a particle can move forward (with probability q) by one 
lattice specing if, and only if, the target site is empty. 
Updating can be implemented either in parallel or in 
random-sequential manner; properties of the model de- 
pends on the updating scheme. For a finite system, either 
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periodic boundary conditions (PBC) or open boundary 
conditions (OBC) can be imposed. When the boundaries 
are open, at every time step, a particle can enter the sys- 
tem (and occupy the site j = 1) with probability a, if 
the site j = 1 is empty. Similarly, under OBC, a par- 
ticle occupying the site j = L can exit with probability 
(3. TASEP under OBC shows an interesting phase dia- 
gram and is the prototype for so-called boundary-induced 
phase transitions Q. To our knowledge, TASEP is the 
simplest model of a system of interacting self-propelled 
particles on a discrete lattice. It has been extended in 
several ways to formulate "particle-hopping" models for 
capturing various interesting aspects of vehicular traffic 

i- 

Our aim is to analyze molecular motor traffic from the 
perspective of vehicular traffic. Therefore, let us first list 
some of the important quantities used in traffic science for 
characterizing traffic flow In the "particle- hopping" 
models of vehicular traffic each vehicle is represented by 
a particle. The number of particles leaving a detector 
site per unit time is defined as the flux and the relation 
between the average flux and the number density of the 
particles is called the fundamental diagram. 

In the particle- hopping models of vehicular traffic 0], 
the time interval in between the departure of successive 
particles from a detector site is defined as the time head- 
way (TH) . Let us start our clock as soon as a particle (let 
us label it by the integer j = 0) just leaves the detector 
whose location is fixed. Suppose, the next n particles 
leave the same detector at times ti,t2, ■■■tn, respectively. 
Then, the corresponding THs are ri — ti, T2 — t2 — ti, 
T3 = ^3 — t2vj Tn = tn — tn-1, respectively. The aver- 
age flux (averaged over the time interval i„) is J = n/i„. 

the mean 



Since YTj=i ^3 = I Ej=i | /n = 1/ J, i.e 
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TH is the inverse of the average flux; therefore, the dis- 
tribution of the THs contains more detailed informations 
on the flow characteristics of traffic than what is available 
from average flux 

The average density-proflle of the particles is yet an- 
other quantitative characteristic of TASEP and TASEP- 
like models. A more detailed characterization of the 
spatial-organization of the particles is possible in terms 
the distribution of the distance-headways (DHs) where 
DH between two successive particles is defined to be the 
number of empty sites in between them [isj . 

In our earlier pap ers on models of molecular motor 
traffic [13, m, [i^IitI. [TsI . [l9j we reported results on aver- 
age flux under both PBC and OBC as well as the average 
density profiles and the phase diagrams under OBC. In 
a few cases [l^, [l^ , we have also reported distributions 
of TH under OBC. In this paper, for each example of 
molecular motor traffic, we briefly summarize the model 
and the fundamental diagram before presenting new re- 
sults on the distributions of THs. We also briefly discuss 
our new observations on the distributions of DHs in these 
models. The TH distributions have important biological 
implications and have been discussed in detail elsewhere 
[la , \W( . Three examples of molecular motor traffic are 
considered in the sections HIl IIIII and IIVI The DH distri- 
butions in the models of molecular traffic are discussed 
briefly in section [V] Our summary and outlook is pre- 
sented in section IVll 



II. 



TRAFFIC OF KINESINS ON 
MICROTUBULE TRACK 



Microtubules and filamentous actin serve as the tracks 
for cytoskeletal motors and both are made of proteins 
[l[. The members of the kinesin and dynein superfam- 
ilies of cytoskeletal motors move on microtubules (MT) 
whereas those of the myosin superfamily move on actin 
filaments. These motors run on chemical fuel in the sense 
that the mechanical energy required for their movement 
is supplied from the energy released by the hydrolysis 
of adenosine triphosphate (ATP); the products of the 
hydrolysis reaction being adenosine diphosphate (ADP) 
and an inorganic phosphate 0|. 

Cytoskeletal motors can attach at any motor-binding 
site on a MT track and can also detach from the track 
[l| . This distinct feature of the cytoskeletal motors was 
captured in the older TASEP-type models of molecular 
motor traffic S [HI [H, US, Hi, [il by adding Langmuir- 
like attachment and detachment processes to the updat- 
ing rules of TASEP. More specifically, in those models, a 
particle is allowed not only to hop forward, but also to 
"attach" to any empty lattice site (with rate uja), and 
"detach" from an occupied site (with rate ujd)- In a 
well-defined regime of parameters, this model exhibits 
co-existence of high-density and a low-density regions 
which are separated by a domain wall (shock); this type 
of spatio-temporal organization is interpreted as a jam 



in molecular motor traffic. But, since most of the effects 
of each mechano-chemical cycle of a motor is captured 
in these models through a single effective hopping rate 
q, it is difficult to make a direct quantitative comparison 
with experimental data which depend on such chemical 
processes. 

We have focussed attention on a family of single- 
headed kinesins, called KIFIA. In order to develope a 
model for cytoskeletal motor traffic, which would not suf- 
fer from the limitations of the earlier TASEP-like models, 
we have incorporated the essential steps of the mechano- 
chemical cycle of individual KIFIA motors in our model. 
Each biochemical cycle of a KIFIA motor consists of a se- 
quence of four states, namely, kinesin (K), kinesin bound 
with ATP (KT), kinesin bound with ADP and phosphate 
(KDP) and, finally, kinesin bound with only ADP (KD). 
The motor binds strongly to the MT track in both the 
states K and KT; the state KDP has very short life time 
and KD binds weakly to the track. Therefore, at each 
spatial location in our simplified model, a KIFIA is al- 
lowed to exist in one of the two distinct "chemical" states 
depending on whether it is bound strongly or weakly to 
the track; these two chemical states are denoted by the 
symbols S and W, respectively. 

The allowed transitions and the corresponding rate 
constants are shown in figdJ The rate constants LUa and 
LU(i account for the attachments and detachments of the 
motors. The rate constant cuh corresponds to the unbi- 
ased one-dimensional Brownian motion of the motor in 
the state where it is weakly bound to the MT track. The 
rate constant ujh is associated with the process driven by 
ATP hydrolysis which causes the transition of the motor 
from the state S to the state W. The rate constants ujf 
and utg, together, capture the Brownian ratchet mecha- 
nism of movement of a KIFIA motor As in the 
earlier TASEP -typ e models of cytoskeletal motor traffic 
[13, mi, m, mil HI, none of the lattice sites is allowed 
to be occupied by more than one motor at a time. 

Let Si {t) and Wi (t) denote the probabilities of finding 
the motor at location i in the strongly- and weakly-bound 
states, respectively. The master equations for Si and Wi 
are given by [l3| 



dSj 
'dt 

<m 

dt 



LOsW^+L0fW,^l[l-S,~Wi) 

+UJa{l- Si~Wi) ~ LOdSi" UJhSi (1) 
UJhS^ + UJb{W^-i + Wr+i){l ~S,- Wi) 

-io,W, - ujfW^il - S,+i - W,+i) 
-L0bW,{2 - S,+i - W,+, - S,-i - iy._i).(2) 



Solving these equations in the steady state under PBC, 
we get the flux 



J = 



LOh + +u;f{l - p) 



(3) 
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FIG. 1: Schematic description of the mechano-chemical cycle 
of a single-headed kinesin motor KIFIA in our model. The 



equispaced sites labelled by the integers 



de- 



note the binding sites of the motor on the microtubule (MT) 
track. The encircled symbols S and W denote the two "chem- 
ical" states of the motor in which it is, respectively, strongly 
and weakly bound to the track. The allowed transitions are 
indicated by the arrows and the symbols accompanying the 
arrows are the corresponding rate constants. 



where 

p = s + w = . (4) 

is the steady-state number-density of the motors on the 
MT track, with K = ujd/oja, = 'jJh/oJf, fls = ojs/cuf, 
and 

D = Afi^Kii + K) + {n,, + ns + {fls - l)Kf. (5) 
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FIG. 2: TH distributions in the model of KIFIA traffic un- 
der OBC. The three curves correspond to three different val- 
ues of uJa (in the units of s~^). The other parameters are 
LOd = 0.1s"\ ujh = 100s"\ LOs = 145s"\ Ljf = 55s~\ 
1.^6 = 1125s-\ a = 50s-^ and /3 = 700s~\ 

In this model, departure of a motor from a site can take 
place in two different ways: either by hopping to the next 



site or by detachment from the track. Therefore, we first 
modify the original definition of TH, which we presented 
in the introductory section, for our model of KIFIA un- 
der OBC. We interpret the attachments at the site i = 1 
and detachments at i = L as consequences solely of for- 
ward hoppings of the motors, rather than manifestations 
of Langmuir-likc kinetics. The time interval in between 
the exit of two successive motors from the site i = L is 
identified as the corresponding TH. 

A few typical TH-distributions under OBC in our 
model of KIFIA traffic are plotted in figH] [26]. A higher 
UJa leads to a higher average density of motors on the 
track which, in turn, reduces the most probable TH as 
long as P remains sufficiently high. This trend of varia- 
tion is consistent with that of the average flux with tOa 
[3]. Moreover, a wider distribution at smaller uJa 
indicates stronger fluctuations in the THs at a lower av- 
erage density. In principle, these new theoretical pre- 
dictions can be tested by carrying out in-vitro experi- 
ments with fluorescently-labelled KIFIA molecules using 
single-molecule imaging techniques. Systematic study of 
the variation of the width of the TH distribution with 
the control parameters of the experiment will provide 
deep insight into the nature of "noise" in intracellular 
cytoskeletal motor transport. 



III. TRAFFIC OF RNAP MOTORS ON DNA 
TRACKS 

Polymerization of a mRNA from the corresponding 
DNA template is carried out by a motor called RNA poly- 
merase (RNAP) and the process is called transcrip- 
tion. To our knowledge, a.11 the models of transcription 
reported earlier [11, S [sil 111, IM H ISi, H llE III 
capture only the stochastic mechano-chemistry of the in- 
dividual RNAP motors. Cooperation and collisions be- 
tween RNAP motors is known to have non-trivial effects 
on the rate of transcription [39, 40, '4lj, Moreover, 
the possibility of the formation of queues in RNAP traffic 
has also been explored (isj. 

The interactions of RNAPs in transcriptional interfer- 
ence [3] is a well known phenomenon and it has also 
been modelled quantitatively [11]. However, instead of 
studying interactions of RNAPs during the transcription 
of different genes, we have modelled the steric interac- 
tions of RNAPs which are simultaneously involved in the 
transcription of the same gene. 

The model is described schematically in figlH Each 
of the lattice sites corresponds to a single nucleotide on 
the DNA template. Each successful addition of a nu- 
cleotide to the elongating RNA is accompanied by a for- 
ward stepping of the RNAP. A single mechano-chemical 
cycle of the RNAP during this elongation stage consists 
of several steps of which the major ones are as follows: (i) 
Nucleoside triphosphate (NTP) binding to the active site 
of the RNAP when the active site is located at the tip of 
the growing RNA transcript, (ii) NTP hydrolysis, (iii) re- 
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FIG. 3: A schematic representation of the mechano-chemical 
cycle of each RNAP in our model. The equispaced sites la- 
belled by the integers + ... denote the nucleotides 
on the template DNA track. No PPi is bound to the RNAP 
in the state 1 whereas the pyrophosphate (PPi)-bound state 
of the RNAP is labelled by the index 2. The allowed transi- 
tions denoted by arrows and the corresponding rate constants 
are also shown. 



lease of pyrophosphate {PPi), one of the products of hy- 
drolysis, and (iv) accompanying forward stepping of the 
RNAP. Since PPi-release is known to be the rate-limiting 
step, we consider only two distinct chemical states fi of 
the RNAP; (1 = 1 refers to the state in which the RNAP 
is not bound to any PPi whereas /i = 2 corresponds to 
the state with bound PPi. 

The processes corresponding to the rate constants uji2 
and W21 are dominated by PP^ -release and its reverse 
reaction. The symbol oj^i is the rate of the successful ad- 
dition of an NTP, catalyzed by the RNAP, whereas U1I2 is 
that of the reverse reaction. The remaining four rate con- 
stants correspond to polymerization/dcpolymerization of 
the RNA, by one monomer, unaided by the RNAP. Since 
premature detachment of a RNAP from its track is a very 
rare event during the elongation of the growing RNA, we 
do not allow such processes in our model. 

The values of some of the rate constants used in our 



numerical studies are as follows: uji2 = 31.4 s 



0.21 s-\ 
0.31 X 10" 



4.66 X IQ- 



9.4 s- 



0.063 s~ . The numerical values 
of the other rate constants are given in the captions of 
the figures [H \E\ and [H 

Unlike the cytoskeletal motors, a single RNAP is so 
large that it can simultaneously cover r successive nu- 
cleotides on the track. In our terminology, a site is occu- 
pied by a RNAP if it coincides with the leftmost of the r 
sites representing that RNAP while the next r — 1 sites 
on its right are said to be covered by the same RNAP. Ir- 
respective of the actual numerical value of r, each RNAP 
can move forward or backward by only one site in each 
time step, if demanded by its own mechano-chemistry, 
provided the target site is not already covered by any 
other RNAP. The total number of RNAPs on the DNA 
template is denoted by the symbol N. Thus, p = N/L is 



the number density of the RNAPs. The coverage density 
is defined by pcov = Nr/L which is the total fraction of 
the nucleotides covered by all the RNAPs together. 
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FIG. 4: Distributions of the time-headways in the model of 
RNAP traffic under PBC. The three curves correspond to 
three different coverage densities of the RNAP motors, all for 



fixed Lu. 



f 



100 and LJ21 = 1.0. 



Let Pf^{i,t) denote the probability that there is a 
RNAP at the spatial position i and in the chemical state 
/i at time t. Let Q{i\j) be the conditional probability 
that, given a RNAP at site i, site j is empty. Note that, 
if site i is given to be occupied by one RNAP, the site i — 1 
can be covered by another RNAP if, and only if, the site 
i — r is also occupied. In the mean-field approximation, 
the master equations for P^(i,t) are given by [16] 



dPiii,t) 
dt 



= 4i Pl(^-l,0 Q(t-l\i-l + r) 

+ Pi(i + l,t) Qjt + l-rU + l) 

+ Lol^ P2{i + l,t) Q{i + l-r\i_±l) 

+ OJ12 P2{i,t) - UJ21 Pi{i,t) 

- (cjfi+w^i) Pl(^,^) Q(i|2 + r) 

- uj'li Pi{t,t) Q{i~r\i) 



(6) 



dP2{i,t) 
dt 



= ^4 P2{i~l,t) Q(i~l\i-l + r) 

+ UJI2 P2{i + I, t) Q{i + I - r\i_±X) 

+ ujj, Pi{i-l,t) Q(i-l\i-l + r) 

+ UJ2lPl{i,t) - UJ12 P2[i,t) 



(7) 
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In the steady state under PBC, 

'uJi2 + ijj\2Q 



Pi 
P2 



where 



f7j + Vl^Q 



fij; = tjJl2 + W2I 
h 



(8) 



21 



^12 



and Q is given by 

Q(i|i + r) = Q{i\ i + r ) = 



1 — pr 



1 + p - pr 

The corresponding steady-state flux is given by 
J = Qi PiQ+ Q2 P2Q 

= {ni Pi+ n2 P2 



(9) 



(10) 



1 - Pc 



1 



Pc 



(11) 



where 



^2 = "^22 ^ ^12 ^ ^22 1 



(12) 



are two effective forward hopping rates from the states 1 
and 2, respectively. 
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FIG. 5: Distributions of the time-headways in the model of 
RNAP traffic under PBC. Different curves correspond to dif- 
ferent values of cu^i, all for the fixed coverage density 0.5 and 
UJ21 = 1.0. 



(i.e., forward hopping) of the successive motors from the 
same site. 

Typical TH distributions of the RNAPs in this model 
under PBC are plotted in figs HI El El for biologically rele- 
vant sets of values of the parameters |43| ; the correspond- 
ing results under OBC have been reported elsewhere p^ . 
In these parameter regimes, because of the "chemical 
transition" between the states 1 and 2, the minimum 
TH is non-zero, i.e., P{t) ~ for all r < T,„i„. The 
magnitude of Tmin as well as the most probable TH are 
practically independent of the coverage density. More- 
over, both Train and the most probable TH decrease with 
increasing Wji, (i.e., with increasing concentration of the 
monomeric subunits of the growing mRNA). In contrast, 
Tmin and the most probable TH increase with increasing 
0^21 which tends to suppress forward movement of each 
RNAP along the main pathway (see figE]). The trends 
of variation of the most probable TH with the model pa- 
rameters are consistent with the corresponding trends of 
variation of the average flux with the same parameters. 

Very recently we have proposed [l^ that the width 
of the TH distribution, as defined above in the context 
of RNAP traffic, can serve as a good quantitative mea- 
sure of noise in transcription. Recent progress in imaging 
techniques has made it possible to monitor the synthesis 
of successive individual RNAs in a living cell ml . U sing 
these techniques it has been demonstrated |47|. [isl liof 
that synthesis of RNAs take place in "bursts" [5c|. The 
time series of the events correspoding to the TH distri- 
butions shown in figs|4l[5]and[6]do not account for such 
"bursts". But "bursting" is observed when our model 
is extended by incorporating the processes of "swtiching 
on" and "switching off" of the gene [sTj . 



o 




2.0 3.0 
Time headway (s) 



5.0 



Since, unlike vehicles in highway traffic, an RNAP mo- 
tor can step backward, TH under PBC must be computed 
carefully as follows [45| : we label the motors sequentially 
in the beginning. Because of the PBC and because of the 
impossibility of overtaking, the sequence of the labels re- 
main unchanged during the time evolution of the system. 
TH is defined as the time interval between the departure 



FIG. 6: Distributions of the time-headways in the model of 
RNAP traffic under PBC. Different curves correspond to dif- 
ferent values of UJ21, all for the fixed coverage density 0.5 and 
ioli = 100. 
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IV. TRAFFIC OF RIBOSOMES ON MRNA 
TRACKS 



Synthesis of a protein from the corresponding mes- 
senger RNA (mRNA) is carried out by a motor called 
ribosome [H^l and the process is referred to as trans- 
lation of genetic code. Usually, many ribosomes move 
simultaneously on a single mRNA strand while each syn- 
thesizes a separate copy of the same protein. Strictly 
speakin g, a ribosome is neither a particle nor a hard rod 
[HS, [13, lE S [13, S [Hi, EH; its mechanical move- 
ment along the mRNA track is coupled to its internal 
mechanochemical processes which also drive the synthe- 
sis of the protein. 

Very recently, we have developed a model that not only 
incorporates the inter-ribosome steric interactions, but 
also captures explicitly the essential steps in the intrari- 
bosome chemomechanical processes. This model for ri- 
bosome traffic is very similar to that for RNAP traffic 
except that (i) each lattice represents a codon (a triplet 
of nucleotides) and (ii) the mechano-chemical of a sin- 
gle ribosome (see figUl) is quite different from that of a 
RNAP. 

Each ribosome consists of two subunits. The smaller 
subunit binds to the mRNA track. But, the actual elon- 
gation of the protein takes place in the larger subunit by 
the addition of the successive monomers, called amino 
acid. The operations of the two subunits are coordinated 
by an adaptor molecule called tRNA. It uses the anti- 
codon at one of its ends to decode the genetic instruc- 
tions stored in each codon of the mRNA which serves also 
as the template for synthesis of a protein. Each tRNA 
carries an amino acid at its other end. Correct codon- 
anticodon matching ensures that the correct amino acid, 
as dictated by the mRNA template, is used by the larger 
subunit for elongating the protein by one monomer. 

The main steps in the mechano-chemical cycle of a ri- 
bosome are the following: 

(i) arrival of the correct tRNA (rate constant uJa), (ii) 
growth of the protein by one monomer because of the for- 
mation of a covalent bond with the newly arrived amino 
acid (rate ujg), and (iii) forward movement of the ribo- 
some by one codon which is associated with the hydrol- 
ysis of a guanosine triphosphate (GTP) molecule (rate 
ujh2)- The rate constants ujhi and are associated 
with two steps during which another GTP molecule is 
hydrolyzed. A more detailed description of the model 
is available in ref. The biologically relevant values 
of the parameters, which have been kept fixed through- 
out our studies of the TH distributions, are as follows: 
u)p = 0.0028s-\ k2 = 2.4s-i and ojg = 2.5s-^. 

Under mean-field approximation, the master equations 
for the probabilities Pf_i {i) are given by 



dt 



uJh2Ph{i-l)QC]zlV-^ + r)+uJpP2{i)- 



-UJaPlii) 

(13) 




FIG. 7; A schematic representation of the biochemical cycle 
of a single ribosome during the elongation stage of transla- 
tion in our model The index below the box labels the 
codon on the mRNA with which the ribosome binds. Each 
circle labelled by an integer index represents a distinct state 
in the mechano-chemical state of a ribosome. The symbols 
accompanied by the arrows define the rate constants for the 
corresponding transitions from one state to another. 
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FIG. 8: TH distributions in the model of ribosome traffic 
under PBC. The three curves correspond to three different 
coverage densities of the ribosomes, all for LJa = 2.5s~^ and 



dP2{i) 
dt 



-uJaPiii)-iL0p + u;hi)P2{i) (14) 



dt 



dPS) 
dt 



UJhlP2ii) - k2P3it) 



= k2P3il)-L^gPi{i) 



(15) 



(16) 



dt 



■ ujgPiii) ~ Uh2P5i^)Q{i\^ + r) (17) 
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FIG. 9: TH distributions in the model of ribosome traffic 
under PBC. Different curves correspond to different values of 
uja, all for coverage density = 0.52 and u^i ~ uJh2 = i-^h = 
Ws-\ 



the parameters of the model of ribosome traffic are very 
similar to those observed in the preceeding section in the 
context of RNAP traffic. Neither Tmin nor the most prob- 
able TH depend on the coverage density. But, both Tmin 
and most probable TH decrease with increasing uja which 
is a measure of the abundance of the amino acid subunits 
of the growing protein chain. A similar trend of varia- 
tion of TH distribution is also observed with increasing 
LUh which is a measure of the rate of "fuel" consumption 
(more precisely, rate of GTP hydrolysis). 

The width of the TH distribution in ribosome traffic 
can serve as a mesure of translational noise [l9l |. just as 
that in RNAP traffic has been considered in section IIIII 
and in ref.[16] as a measure of transcriptional noise. A 
comparison between the trend of variation of this noise 
with our model parameters and the corresponding recent 
experimental observations has been reported elsewhere 
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V. DH DISTRIBUTIONS IN MOLECULAR 
MOTOR TRAFFIC 

Because of the possibilities of attachments and detach- 
ments of the motors, the DH distribution is not a suitable 
quantity for characterizing KIFIA motor traffic. There- 
fore, in this section, we shall consider exclusively only 
traffic of RNAPs and ribosomes under PBC; the corre- 
sponding results under OBC are discussed in ref.ji^ and 
ref.[2^. 

The DH distribution in the TASEP with hard rods, 
each of length r, was calculated analytically by Shaw et 
al. [56(1 . Suppose, m denotes the DH. The DH-distribution 
P{m) in this model under PBC is given by [56| 



FIG. 10: TH distributions in the model of ribosome traffic 
under PBC. Different curves correspond to different values of 
uih, all for coverage density — 0.52 and Ua = 2.5s~^. 



In the steady state under PBC, the flux of the ribo- 
somes is given by 



J = 



i^h2p{'^ - pr) 



where 



with 



{1 + p - pr) + nh2{l - pr) 



^h2 = ^h2/keff- 



UJa 



1 

k'2 



1 



1 



(18) 



(19) 



(20) 



TH distributions of the ribosomes in this model under 
PBC are plotted in figs El [H [10] for different sets of pa- 
rameters [26| ; the corresponding results under OBC have 
been reported elsewhere [l^. The qualitative features of 
the TH distributions and their trend of variation with 



(21) 



where ph = 1 — Pcov is the fraction of the system covered 
by the holes, and ps = p+Ph- In the limit r = 1, pcov ~ P, 
Ps ~ \ and, hence, P(m) reduces to the limiting form 
P{m) = p{l — p)"^ which is the well-known mean-field 
estimate of the DH distribution for TASEP [12, d- 

In figlTT]we have plotted the DH distributions in the 
models of RNAP traffic and ribosome traffic under PBC 
where r = 10 for both the models. Interestingly, for any 
given coverage density, the DH distributions in both these 
models follow equation (|21|1 in spite of the differences 
in the mechano-chemical cycle of the individual motors. 
Thus, the DH distributions in these models of molecu- 
lar motor traffic are determined solely by the geometric 
parameters. 



VI. SUMMARY AND CONCLUSION 

In this paper we have modelled molecular motor trans- 
port on filamentary tracks from the perspective of physi- 
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FIG. 11: DH distribution in RNAP traffic and ribosome traf- 
fic. 



equilibrium systems of interacting self-driven particles. 
Moreover, many extensions of TASEP have been used 
earlier to capture various phenomena observed in vehic- 
ular traffic. The interesting quantities which are used for 
characterizing traffic flow include fundamental diagram, 
distributions of distance- and time-headways. In this pa- 
per, for all the models of molecular motor traffic, we have 
not only summarized the earlier results on the funda- 
mental diagrams, but also presented many new results 
on the distributions of DH and TH. The TH distribution 
in traffic of RNAPs and ribosomes yield novel measures 
of intrinsic noise in transcription and translation, respec- 
tively [l^, [l^. Our modeling provide new insight into 
intracellular motor-driven processes by looking at these 
from a novel perspective. 



cists and traffic scientists. These models are essen- 
tially biologically motivated extensions of TASEP. To 
our knowledge, TASEP is the simplest model of non- 
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